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WAVELET TRANSFORM AND DIFFUSION EQUATIONS: 
APPLICATIONS TO THE PROCESSING OF THE 
"CASSINI" SPACECRAFT OBSERVATIONS 



We show that continuous transform with the complex Morlet wavelet is easily per- 
formed if we replace the integration of the fast-oscillation function by the solution 
of the diffusion differential equations. The most important advantage of this ap- 
proach is that the initial data can be represented by non-uniform sample of an 
arbitrary node number. We apply the proposed method to the processing of the 
image of the Saturn A-ring obtained from the Cassini spacecraft. Also we have got 
via the wavelet transform using PDE, the local correlation coefficient of the signal 
and the harmonic function. 



1. Introduction 

As is known, during last 20 years the wavelet transform is one of the con- 
temporary methods for the signal processing. The main advantage of this 
approach is the high localization of the base functions in both, spatial 
and frequent regions (see, e.g., 1 , 2 ). This allows us to analyze effectively 
non-stationary signals and signals with certain singularities. For our con- 
venience, the continuous wavelet transform 



E.B. POSTNIKOV 



Kursk State University, Theoretical Physics Department 
305000, Radishcheva st., 33, Kursk, Russia 
E-mail: postnicov@mail.ru 



A.LOSKUTOV 



M. V. Lomonosov Moscow State University, Physics Faculty, 
119899, Leninskye Gory, GSP-2, Moscow, Russia 
E-mail: loskutov@polly.phys.msu.ru 




(1) 



— OO 



February 2, 2008 16:29 Proceedings Trim Size: 9in x 6in 



Losk' article 



2 



(where the asterisk means the complex conjugation) with the amplitude 
norm 



DC 



t-b 



dt , . 

— = const. (2) 



is the most convenience. 

The wavelet basis which is often used for the signal processing is the 
complex Morlet basis. Its exact form that satisfy the admissible condition 



/oo 
-oo 



has the following form 2 : 



^) = ^=(e-^- e -^)e-V. ( 3 ) 



In this case, the corresponding wavelet-transform w(a, b) plays the role of 
the local spectrum distribution by period a in the neighborhood of the point 
b. This feature allows, on the basis of the instant frequency notion 3 , to 
develop the method of the analysis of the solution for Hamiltonian systems. 
Another approach developed on the basis of the mean-square deviation of 
the investigated series from the model data (but without of the wavelets) 
has been proposed in 4 . In the present paper, we show that these methods 
are similar to each other. 

It should be however noted, that the standard realization of the con- 
tinuous wavelet transform based on the transition to the frequency domain 
and fast Fourier transform (FFT) has certain disadvantages. They follow 
directly from FFT: the initial data should be presented by the sample 2 N 
equispaced nodes. Neglect of this condition leads to the algorithm com- 
plication and the lost of the accuracy. But the data sets obtained via 
dynamical systems and some astronomical observations can be presented 
by sufficiently non-cquispaccd sample. 

Thus, we should rely upon an alternative method of the evaluation of 
the wavelet transform with the real-valued Gaussian family: 

1 d n _ « 2 

^ (0 = ^dr e " T ' 



This approach presents the wavelet-transform as a derivative of the diffu- 
sively smooth initial function 5 : 
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This function can be found as a solution of the diffusion equation with the 
iterative or finite-difference methods. In this case, there are no restrictions 
to the initial date representation. In addition, this approach can be adapted 
for the multidimensional systems. For example, in 6 it successfully used for 
the analysis of 3D turbulent solutions. 

The main purpose of the present paper is to reduce the continuous inte- 
gral wavelet transform with complex wavelet to the solution of the partial 
differential equations (a Cauchy problem) . It allows us to use the developed 
algorithm in the local frequency analysis. This approach (which cannot be 
realized by the standard FFT) additionally makes it possible almost arbi- 
trary scale decompositions. This problem is very important now in applica- 
tion to small-scale structure of the planetary rings, in particular, the Saturn 
rings (see, e.g. a review 7 ). For the first time, analysis the Voyager-2 data 
by the wavelet transform has been undertaken in 8 . The authors using the 
real wavelets, detect the structure of the rings in the noisy images. At the 
same time, the problem related to the local periodicity is still open, but it 
can be resolved by the complex wavelet transform. This is due to the fact 
that in the present the necessary data are available by the Huygens-Cassini 
mission. 



2. Continuous transform with the Morlet wavelet 

For applications, when loq is a large enough, in (3) one can neglect the 
second summand. Thus, in the present Part we will used the following 
simplification: 



V(0 = ^=e--°€ e -V. ( 4 ) 



This expression corresponds to the normalization (2) with const = 
exp(ojQ/2). Let us show that such a choice (owing to the violation of ad- 
missible conditions) makes it possible to use the transformed function as 
an initial condition for the equation 

d 2 d d\ 
a 7^I ~ q !W077 P(0,6) = (5) 



db 2 da db 

which is satisfied by the wavelet-image (I) with the Morlet wavelet (see 9 ). 
Let us write the wavelet-image as a sum of the real and the imaginary parts 



w(a, b) = u{a, b) + iv(a, b). 
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With respect to these variables Eqs.(5) may be presented by the following 
system: 

du d 2 u dv 

d^= a dV +UJ °db (6) 
dv d 2 v du 

To find the corresponding initial conditions let us rewrite the continuous 
transform (1) with the kernel (4) in the form: 



w(a,b)= J f(t f 2 J 2 ^ 2 ' dt. (8) 



As is known, the integral (8) does not depends on the imagine subtrahend 
in the exponent. Also, in the limit a —* the transformation kernel is a 
delta-function. Therefore, we can get initial conditions 

u(0,b) = Re(f(b)), 
w(0,6)=/m(/(6)) 

for the system of differential equations (6)-(7). 

Let us apply the described algorithm to the analysis of radial distribu- 
tions of the matter density in the Saturn A-ring. Here we use the data 
obtained in July 2004 from the Cassini orbiter (image PIA 06091 from the 
collection NASA/JPL/Space Science Institute). This image consists of a 
stripe with the width of 29 pixels in the radial direction. The natural length 
of the image is about 248 kilometers, and the width is about 7.2 kilometers. 
In such an approximation we may neglect the curvature of the ring. The 
function f(t) obtained by the sample average is shown in Fig. la. 

As the signal is rcalvalued, we used the following pair of initial condi- 
tions: u(0, b) = f(b) and v(0, b) = 0. Owing to the bound size, the Cauchy 
problem for the Eqs. (6)-(7) can be replaced by the boundary problem 
of the first kind. For the real component of the function we used the ini- 
tial signal value. For the imaginary component zero boundary conditions 
has been taken. Because it is necessary to get the most good resolution 
for the present numeric sample, the based frequency has been chosen as 
ujq = ir. By the reason that the inverse wavelet transform is not necessary, 
we neglect by the deviation of the wavelet for the admissible conditions. 

The results of numerical analysis, the modulus and the phase of the 
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Figure 1. The radial distribution of particle, modulus and the phase of the wavelet- 
transform. 

wavelet-transform determined as 

\w(a, b)\ — \/ u 2 (a, b) + v 2 (a, b), 6(a, b) = arctan U ^f" 1 ^ , 

v(a,b) 

are shown in Fig. lb, c. In this Figure, the larger values correspond more 
dark regions. 

Let us consider the modulus of the wavelet-transform (Fig. lb) in detail. 
One can easily see that there arc regions with monotonically changing in- 
stant periods. These regions belong to the intervals 0.05 — 0.2 and 0.5 — 0.85. 
In the non-regular region of the ring A periodic components do not revealed. 
Note that to avoid the boundary effect, we do not consider intervals closely 
related to the ends of the sample. 

As is known, values and dynamics of the instant period are determined 
by the modulus maxima lines of the wavelet-transform. Let us extract them 
from the region with the regular image (see Fig. lb). The obtained lines are 
shown in Fig. 2a, b. The form of these lines confirm the understanding of 
the fine stricture as the regions corresponding to the high-order resonances 
with Saturn's moons. In particular, the analyzed image PIA 06091 from 
the collection NASA/JPL/Space Science Institute can be construed as the 
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region between the resonances 11 : 12 Prometheus and 3 : 5 Mimas. Using 
the third Kepler law we may obtain the distribution of the major semi- 
axes corresponding to the Farey sequence from this interval (Fig. 2c, d). In 
this Figure, we clearly see the near resemblance of the found parts of the 
sequence with the line of the instant radial periods. 




Figure 2. Instant periods of the radial particle distribution in the Saturn A-ring (a,b) 
and the Farey sequence (c,d). 



3. The wavelet transform with the Morlet basis and the 
correlation coefficient 

The rigorous definition of the Morlet basis (3) makes it possible to find 
quantitatively the correlation of the instant period with period of the 
monochromatic harmonics. Introduce a new denotation a = a/u;o and 
rewrite the corresponding integral convolution: 



f(t)e 



-.dt. 



(9) 
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At f(t) — 1 the integral (9) has a closed form solution e'»e 2° . Thus, the 
following representation is admissible: 

+ OC 



e 2 — e 



(a -l )t J — 0^ 



V /27TQ! 2 Wq 



(10) 



Therefore, the wavelet transform (1) with the Morlet basis (3) can be rewrit- 
ten as: 



w(a, b) = e 



f(t)e^^ -, 



:dt- 



+ OC - 2 +°° 

/ e^^2=dt I fit) 



Let us compare here the expression in the square brackets with the corre- 
lation coefficient of two random variables x and y: 

R = (xy)- (x)(y) 

where (•) means the averaging operation and = (£ 2 ) — (£) 2 is a dispersion. 
We see that this wavelet transform is equivalent to 

w(a,b)=e-^ [(/(^e^" 1 )*)-^^))^" 1 )*) 

and t are normally distributed sample nodes. This distribution is centered 
in b. 

Let us introduce dispersions: 

2 



+ 00 



(t-b) 2 



ITCt^UJ, 



2„,2 



.dt 



f(t) 



zdt 



Dexp — 



(*-*>): 

2q 2 w 2 

= 2 '(«- 1 )*^ -.dt 



v / 2ira 2 L 



I +oc 

- / «**-'' 



(t-fc) 2 



^27ra 2 c 



The second expression has the following simple analytical solution: 



n — e -^o _ p - 2w o 

^exp — c c 
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Now we are able to find quantitatively the local correlation coefficient of the 
signal /(f) and the function exp(i(a _1 )f) in a neighborhood of the given 
point f = b: 

R a (b) = ^L. (11) 

yf JJfJJ exp 

All the integrals contained the analyzed function are the solution of the 
Cauchy problem for the simplest diffusion equation: 

du{r,b) 2 d 2 u(T,b) 
dr db 2 

These solutions arc taken in time t = uJq/2. As an initial condition we used 
/(f) exp(— (a _1 )f), /(f) or /(f) 2 . To solve this equation there are robust 
methods (see, e.g. 10 ) even on the non-uniform grid. 



4. Conclusions 

Thus, we showed that continuous transform with the complex Morlet 
wavelet is easily performed if we replace the integration of the fast- 
oscillation function by the solution of the differential equations. In spite 
of fact that this method requires more computer resources, it has certain 
advantages. The most important of them is that the initial data can be 
represented by non-uniform sample of an arbitrary number nodes. In ad- 
dition, this algorithm is realized via the standard MATLAB. During the 
calculation, the time variable in the diffusion equation is associated with 
the scale variable of the wavelet transform. Therefore, the choice of a small 
enough step allows us to trace the evolution of the instant period in detail. 

In the present article we apply the proposed method to the processing 
of the image of the Saturn A-ring PI A 06091. This analysis confirm the 
resonance origin in its structure. 

More detail quantitative description is the subject of our following anal- 
ysis. This analysis is possible due to the video material obtained from the 
Cassini spacecraft last half 2004 year. These quantitative estimations can 
be based on results descried in the second part of the article, where we got 
the dependence of the local correlation coefficient of the signal and the har- 
monic function in neighborhoods of the given point. Also, the computation 
of this coefficient may be realized via the solution of diffusion equations. 

The proposed algorithm may be additionally considered as a certain ex- 
tension of the local time-series analysis. Namely, in , a set of sequence with 
the compact support has been advanced as test values. These sequences 
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are obtained by the periodization of sample's parts. In our analysis, as a 
test we propose to use the harmonics. 
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